Breakup of a Stoner model for the 2D ferromagnetic quantum critical point 
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Re-interpretation of the results by [A. V. Chubukov et. al, Phys. Rev. Lett. 90, 077002 (2003)] 
leads to the conclusion that ferromagnetic quantum critical point (FQCP) cannot be described by a 
Stoner model because of a strong interplay between the paramagnetic fluctuations and the Cooper 
channel, at least in two dimensions. 
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Recent experimental observations of superconductivity 
in a close proximity to the ferromagnetic quantum criti- 
cal point (QCP) in a few itinerant ferromagnets such as 
heavy fermion compounds UGe^i, URhGe^ and also in 
ZrZii2' ! have revived an interest to whether the magnetic 
fluctuations in the vicinity of ferromagnetic QCP provide 
a fundamental mechanism for the "p-wave" superconduc- 
tivity. The idea of the exchange by paramagnons turned 
out to be crucial for understanding the mechanism un- 
derlying the superfiuid state in 3 He at low temperatures . 
In a broader context, relevance of paramagnetic fluctua- 
tions to the appearance of superconductivity at the bor- 
der of a magnetic transition has been discussed in— ano— 
where reader will find the history of the problem and the 
comprehensive list of references. 

Starting already with the one of the first papers on 
the subject there were numerous efforts to evaluate 
the superconducting transition temperature mediated by 
strong paramagnetic fluctuations for both ferromagnetic 
and paramagnetic phases&Siiil. Most of these results bear 
the numerical character, and this obscures the fact that 
there is, as we believe, some flaw in the very model. Most 
of the authors deal with the Stoner model where an in- 
stability of the paramagnetic itinerant state comes about 
with the increase of the on-site Hubbard interaction, U, 
that leads to the appearance of ferromagnetism. Varia- 
tion in the value of the Stoner factor governs then changes 
in a value of the Curie temperature by varying an exter- 
nal parameter (in a more general form these ideas can 
be formulated in terms of the Fermi liquid theoryii by 
involving the Pomeranchuk instabilitjiiSii 3 .) . Below we 
will try to demonstrate that such a model is not self- 
consistent, at least for the 2D systems. 

The first question one faces while discussing any phase 
transition is how close one can approach the line of a 
transition. The problem of a singularity, or magnitude of 
fluctuations, looks simpler near the ferromagnetic QCP 
(the ending point of the phase diagram at T = 0) . The 
outcome of the analysis done by Hertz 1 and recently by 
Millie 15 is that near such QCP fluctuations retain their 
mean field character due to an increase in effective dimen- 



sionality to account for an involvement of the frequency 
variable at zero temperature. We argue, that the hy- 
pothesis of the ferromagnetic QCP with the paramagnon 
propagatoi*^ would lead to the developing of supercon- 
ducting fluctuations at such a scale which breaks the va- 
lidity of Hertz analysis far away from the vicinity of the 
imaginable QCP. In other words, a Stoner like ferromag- 
netic QCP is not self-consistent namely because it seem- 
ingly leads to such strong pairing fluctuations that make 
incorrect independent analysis of the spin "zero sound" 
and Cooper channels. Below, we would like to prove our 
point by re-interpreting the results of the recent paper 
by Chubukov et aiw&. 

The emphasis in Ref— has been put to demonstrate 
that superconducting transition near the ferromagnetic 
QCP may turn out to be of the first-order. The authors 
oflfi used the standard ansatz for the longitudinal fluc- 
tuations propagator—. To get rid of the so-called non- 
adiabatic corrections and to reduce the gap equations to 
the well known form of the strong coupling ones (" Eliash- 
berg" equations), the authors made an assumption that 
the interaction of electrons with the spin fluctuations is 
weak. A minor change below makes it more convenient 
to overview the physical picture of their model as a whole 
without resorting to numerical calculations. 

Let us introduce the exchange part of interaction be- 
tween the two electrons, I(r± — r 2 ) as: 



I(ri - r 2 ) = -I(ri - r 2 ) (<ri • <t 2 ) 



(1) 



and assume that I{r-y — r 2 ) has the ferromagnetic sign 
and bears a long-range character. After summing up all 
diagrams in the zero-sound channel the same way it was 
done ini4, one arrives to the following longitudinal spin- 
spin fluctuation propagator: 



I(q) 



1-2 v{e F )I{q) 



1-a 



(2) 
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where I{q) is the Fourier component for the interaction 
QJ, w n = nT(2n + 1) is the Matsubara frequency, v{ep) 
is the density of states at the Fermi level and the value of 
a depends on the dimensionality of the system: a = ir/2 
for 3D and a = 1 for 2D. 

Cancellation in the denominator of J5J at q — ► 
and Lu n /vpq <C 1 leads to the Stoner criterion: 1 — 
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2 i/(s F )I(0) — > 0. The factor in the square brackets in de- 
nominator in @ is nothing but the electron polarization 
operator II(q; uo n ): 

r%; u n ) = £ f d d l G(q + l;u n + fi,) • G(l; n,). 

proportional to the generalized electron spin suscepti- 
bility and calculated at small enough q in the limit 
<C v fQ ( a t this relation w n and g appear in all equa- 
tions below). In the opposite limit Tl(q;ui n ) decreases to 
zero, and the Stoner like enhancement rapidly disappears 

at U) n ~V F q. 

According toi^, Eq. © for the ferromagnetic fluctua- 
tions does not experience renormalization at T = and 
1 - 2 v(e f )I(0) -> 0. To agree the form of Eq. © with 
the similar expressions ini& we assume that I(q) rapidly 
decreases with an increase in q. To be more specific, we 
accept the following notation: 

2I{q)u{e F ) = I a cxp [-q 2 /(p F $) 2 ] , (3) 



where 

t? = {RvfV 1 < 1 (4) 

is the parameter sharing the long range character of Q , 
which in turn leads the small angle scattering to prevail 
in Eq. @. The interaction is isotropic in the spin 
space. It is straightforward to account in Eqs. (f 1121) for 
the presence of the magnetic anisotropyS. 

As it was first pointed out ini£, spin fluctuations pro- 
duce two effects in case of the phonon-mediated (s - wave) 
pairing: they add to renormalization of the electron self- 
energy and they provide the pair breaking mechanism 
for the S = Cooper pair. Pair breaking effects are 
basically the same even for a triplet (S = I) pairing, 
however, the exchange between two electrons by longi- 
tudinal paramagnetic fluctuations leads to the attrac- 
tive interaction in the triplet (S — 1) channel. For 
the electron's Green function in the paramagnetic phase, 
G^ 1 (p; uj n ) = iu„ - e(p) - T,(p;u n ), we write: 



S(P;W„) = J^y T Yl j dd Pl I (P - Pl^n~ ^n')G(p 1 ;LJ n/ ), 



(5) 



where d is the dimensionality of the problem, and the coefficient c in © depends on whether the exchange is isotropic 
(c = 3) or a strong magnetic anisotropy is present (c = 2 for "easy plane" and c = 1 for "easy axis"). For the sake 
of simplicity we take c = 1 ("easy axis") in order to avoid additional complications related to the possibility of a 
first-order superconducting transition. 

The superconducting order parameter A a p(p; oj n ) is chosen below in the form: 



A a/3 (p;w„) = i[(a • d(p; u„))a y ] a0 

and near transition the vector d(p;uj n ) satisfies the linear equation: 

T 



d(p;uj n ) = 



(2ir) d 



d PiI(P — Pi;w„ - u n ') G(pi;ov) G(-px;-u n i)d(pi;u n i) 
I 



(6) 



(7) 



which would serve to determine the dependence of the 
superconducting transition temperature on the proximity 
to the QCP in Eq. @. According toi&, the vicinity of 
the ferromagnetic QCP where onset of superconductivity 
acquires a non-perturbative character is much broader 
in 2D then in 3D. Consequently, our discussion will be 
restricted to the 2D situation only. 

Let us express vector q in Eq. through the angle 
of scattering along the Fermi surface, ^)« 1: 



2 2 
PF V 



(8) 



In Eqs. (|5I7|) one can integrate over the energy variable 
component, £i = v F (jp\ —Pf)- Dependence on | p | in the 
order parameter A a f}(p;uj n ) can be neglected because of 



smallness of the scattering angle, (p: 

d{p;u n ) => d(u) n ). (9) 
Introducing the notation for proximity to the QCP 



1 - I Q = t < 1 
after simple calculation Eq. J7J acquires the form: 



(10) 



where 



< J >" = 27 



iu n > - E(uv) 



ds 



^d(uv), (11) 



(12) 
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comes about after integrating J5J) over ip and making use 
of Eqs. (l.'il ll-SIl in the expansion: 



\-2I{q)v(e F ) ~r+(p/0)' 



(13) 



In equations above we have neglected all the vertex 
corrections to the bare vertices. With the help of the 
explicit expressions for I(p;u> n ), it will be easy to ver- 
ify that the higher order corrections are small provided 
that our parameter i? above is small. Straightforward cal- 
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FIG. 1: Vertex correction to the interaction. Dashed and 
solid lines are the paramagnon and fermionic propagators, 
respectively. 

culation of the first vertex correction to the self-energy 
E(p;cj„) (J5J shown in Fig. ^ 

£p(i) _ [sign(uj m + fij) - sign(w m )] x 



dip 



gives the following order of magnitude estimate: 

<5r« ~ (d/V^) ■ (sii/v F k) 

In notations (|12|1 denominator in ljll|l is: 

iui n - £(u%) = iuj n + iitT ^ (-0"v 



(13') 



(14) 



a; / — — a;,. 



Substitution s = z^fr transforms l|12(l into: 

oo 

.—. ( d \ f zdz 



with 



(15) 



(16) 



cj = e F $ 4 /32. 
Frequency dependence in (|T2|l serves as the cut-off: 

(7) w « l/3x/3> /^) 1/3 (17) 

for 



w > wo 



(17') 
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FIG. 2: Numerical solution of Eq. JTTJ for T c as a function 
of t for the various values of magnetic anisotropy constant, 
c (see inset). Maximum in T c for easy plane and isotropic 
cases corresponds to the onset of the first-order transition 
into superconducting state. 

Let us consider the paramagnetic phase (r> 0) far 
enough from the QCP. Renormalization of the self-energy 
part T,(oj n ) is small, so that Eq. (|llfl converts into the 
weak-coupling problem: 



1 = 



/ 9 



(18) 



where ui is of the order of u> (^r^J ( see Eq. (|17'|) h 

When y/r is further decreased, so it becomes comparable 
to the strong-coupling regime sets in with the cut- 
off frequency uj ~ u>q in i|18|) . This leads to the new 
and important changes in £(w). As it has already been 
discussed mi&, one has the following asymptotic behavior 
of the self-energy part: 



E(w) = 



I 3\/3 

V 4 



2/3 



uj < 2e F i9r 3 / 2 , 
uj > 2e_F^r 3 / 2 



(19) 

The second asymptotic would signify appearance of the 
non-Fermi liquid regime in the close proximity to the 
QCP. It is seen, for fixed that the overall behavior 
of the critical temperature is defined by the value of the 
"coupling" constant, A: 



(20) 



As it is readily seen from Eq. . T c ~ ujq at A ~ 1. 

Equation l|ll|) for T c as a function of r was solved nu- 
merically with the solution shown on Fig. At r = 
we obtain: 



T c ~ 2.5w - 0.08£Ftf 4 . 



(21) 



Thus we conclude that T c is finite at r = and reaches 
the energy scale of the order of ui already at A ~ 1 . The 
latter is also true for other cases shown on Fig. 2. 
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Let us therefore keep A ~ 1 and start gradually in- 
creasing value of the model parameter At d — > 1 we 
should return back to the Stoner model with ordinary 
local interaction which is discussed ini4. We see from 
(|13'|) , that non-adiabatic corrections remain of the order 
of one at A ~ 1, $ ~ 1 and, hence, one expects that their 
exact treatment would not change qualitatively the es- 
timate T c ~ Ep from Eq. 1|21|) . On the other hand, at 
A ~ 1 and -Q — > 1, r = i? 2 /16A 2 ~ 1, and one cannot come 
close to QCP without forming a superconducting ground 
state, which in turn would change the polarization opera- 
tor |J2j. Obtaining such a high values for T c energetically 
so far away from the originally accepted position of the 
QCP (t = 0) shows the intrinsic contradiction of the local 
Stoner model^. At such a high energy scale the assumed 
proximity to a " QCP" seems to be irrelevant for physical 
properties of the system. Mathematically, large T c ~ ef 
or/and large A lead to the change in the polarization 
operator Eq. JJJ to account for an interplay between 
the zero-sound and superconducting channels. The po- 
larization operator J2]) is modified due to the presence 
of anomalous Gor'kov functions F, in superconduct- 
ing state at T = 0. For ■d ~ 1 and A ~ T c ~ e F this 
introduces such a change in polarization operator that 
significantly reduces cancellation in the Stoner factor. 

Remember now that at u> — and vf<1 — ► expression 
used to describe a behavior of magnetic susceptibil- 
ity near a ferromagnetic QCP. Similarly, modified polar- 
ization operator is proportional to electronic susceptibil- 
ity in superconducting state. The latter would not go to 
zero at T = as it does for the s-wave pairing, even 
though it does not equal to normal susceptibility neither 
in any triplet state. The Stoner cancellation does not 
occur. Therefore, once superconductivity H18|) arises at 
t > in the framework of the model with d <C 1, at 
T = the ground superconducting state as a function 
of external parameter continues to be stable while en- 
tering into the ferromagnetic state (r < 0) well beyond 
"QCP". This suggests a first order like phase compe- 
tition between superconductivity and ferromagnetism at 
low temperatures. 

All that has been said above poses a few questions. 
First, numerical calculations (see Ref&2iiii) for T c mak- 



ing use of an exchange by longitudinal spin fluctuations 
give rather low values of critical temperature compared 
to the values of the bandwidth or the Fermi energy with- 
out special assumption of small angle scattering (the au- 
thors have been solving basically the same equations, i.e. 
no vertex corrections have been included). We believe, 
this is a result of some numerical smallness, such as the 
factor 1/32 in Eq. Q16p. This smallness may restore 
QCP within some vicinity. Indeed, the low values of T c 
has been experimentally observed in a number of systems 
among which are Sr2Ru04, URhGe or 3 He. In addition, 
in 3D the strong coupling regime of the above model sets 
on only in very close proximity to the QCP—. Note that 
the model itself does not determine the spatial symme- 
try of the triplet Copper pair wave-function: exchange 
by paramagnons is not the only interaction between the 
electrons in the Cooper channel. Thus, the paramagnon 
contribution in our model comes together with other in- 
teractions, r; ~ 1: 

(3D): ~$ 2 \ln(T)\+T i; (2D): ~^+T h (22) 

where Ti may be positive (here T; is a proper harmonics 
designed by the exact symmetry of the superconducting 
pairing). Onset of attraction in ill'L'l) takes place closer to 
t = and an effective interaction 122|) remains reason- 
ably weak (~ 1) in its vicinity^. 

To summarize, we have shown that competition 
with the superconductivity channel makes the Stoner 
model for ferromagnetic quantum critical point not self- 
consistent. FQCP can be realized only due to the pres- 
ence of a numerically small parameter or other repulsive 
interactions in the triplet channel that weaken the at- 
traction mediated by paramagnons, at least in 2D. 

The authors thank A. Chubukov, A. Finkel'stein and 
D. Morr for fruitful discussions. 

This work was supported by DARPA through the 
Naval Research Laboratory Grant No. N00173-00-1-6005 
(M.D.) and by the NHMFL through the NSF coopera- 
tive agreement DMR-9527035 and the State of Florida 
(L.P.G.). 



1 S. S. Saxena et al, Nature 394, 39 (1998). 

2 D. Aoki et al., Nature 413, 61 (2001). 

3 C. Pfieiderer et al, Nature 412, 58 (2001). 

4 A. J. Leggett, Rev. Mod. Phys. 47, 331 (1975). 

5 N. D. Mathur et al, Nature 394, 39 (1998). 

6 A. Chubukov, D. Pines and J. Schmalian, in The Physics 
of Conventional and Unconventional Superconductors, ed. 
by K.H. Bennemann and J.B. Ketterson (Springer- Verlag, 
2002) 

7 D. Fay and J. Appel, Phys. Rev. B 22, 3173 (1980). 

8 Z. Wang, W. Mao, and K. Bedell, Phys. Rev. Lett. 87, 
257001 (2001); K. B. Blagoev, J. R. Engelbrecht, and K. 



Bedell, Phys. Rev. Lett. 82, 133 (1999). 
9 R. Roussev and A. J. Millis, Phys. Rev. B 63, 140504 
(2001). 

Ph. Monthoux and G. G. Lonzarich, Phys. Rev. B 59, 
14598 (1999). 

1 I. E. Dzyaloshinskii and P. S. Kondratenko, Sov. Phys. 
JETP 43, 1036 (1976). 

2 I. Ia. Pomeranchuk, Sov. Phys. JETP Lett. 8, 361 (1959). 

3 A. A. Abrikosov and I. E. Dzyaloshinskii, Sov. Phys. JETP 
35, 535 (1959). 

4 J. Hertz, Phys. Rev. B 14, 1165 (1976). 

5 A. J. Millis, Phys. Rev. B 48, 7183 (1993). 



5 



16 A. V. Chubukov, A. M. Finkel'stein, R. Haslinger, and D. 
Morr, Phys. Rev. Lett. 90, 077002 (2003). 

17 N. F. Berk and J. R. Schrieffer, Phys. Rev. Lett. 17, 433 
(1966). 

18 adding T; > breaks the logarithmic self-consistance due 



to the form of the kernel 11711 and the self-energy part 
E(u) dJl. Non-zero Fi in (1221 narrows the vicinity of QCP 
where solutions for T c shown on Fig. 2 apply. 



